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Background and purpose — The continuously increasing demand 
for joint replacement surgery in the past decades imposes higher 
constraints on the budgets of hospitals and healthcare providers. 
We undertook an analysis of historical trends in total hip replace- 
ment performed in Sweden between 1968 and 2012 in order to 
provide projections of future demand. 

Data and methods — We obtained data on total hip replace- 
ments registered every year and on the evolution of the Swedish 
population between 1968 and 2012. We assumed the existence of a 
maximum incidence. So we adopted a regression framework that 
assumes the existence of an upper limit of total hip replacement 
incidence. 

Results — We found that the incidence of total hip replacement 
will continue to increase until a projected upper incidence level of 
about 400 total hip replacements per 10 s Swedish residents aged 
40 years and older will be reached around the year 2107. In 2020, 
the estimated incidence of total hip replacement will be 341 (95% 
prediction interval (PI): 302-375) and in 2030 it will be 358 (PI: 
317-396). Using official forecasted population growth data, about 
18,000 operations would be expected to be performed in 2020 and 
20,000 would be expected to be performed in 2030. 

Interpretation — Growing incidence, population growth, 
and increasing life expectancy will probably result in increased 
demand for hip replacement surgery. Our findings could serve as 
a basis for decision making. 



Osteoarthritis is a leading cause of disability in the elderly. 
Patients with osteoarthritis contribute to higher annual health- 
care-related costs than those without (Kotlarz et al. 2009, 
Berger et al. 2011). Increase in life expectancy and ageing 
populations are expected to make osteoarthritis the fourth 
leading cause of disability by the year 2020 (Woolf and 
Pfleger 2003). Thus, the demand for joint replacement surgery 



due to osteoarthritis will grow, imposing higher constraints on 
the budgets of hospitals and healthcare providers. Decision- 
makers and insurers require reliable and accurate projections 
of future demand, for proper budgeting and fair allocation of 
resources. Increased demand for joint replacements also cre- 
ates an increased need for surgeons trained and experienced in 
joint replacement surgery, which in turn requires educational 
planning at universities and hospitals. 

Recent studies have forecast a continuously increas- 
ing demand for both hip and knee replacement (Kurtz et al. 
2007a). Primary osteoarthritis, the leading determinant of 
total hip replacement (THR), has a more or less age-specific 
occurrence and a relatively low incidence in the general popu- 
lation. Thus, the existence of a maximum incidence is a natu- 
ral assumption. This must be built into any analysis for the 
statistical model to be adequate (Turner et al. 1961). Analysis 
of incidence rather than the total number of THRs allows not 
only projections of future demand, but also an estimate of the 
highest expected incidence of total hip replacement. In this 
study, we undertook an analysis of historical trends in THR 
performed in Sweden between 1968 and 2012. Our first aim 
was to forecast the incidence of THR operations in Sweden 
in the coming decades. The second aim was to estimate the 
expected maximum incidence of THR operations in Sweden. 



Data and methods 

Ethics 

Ethics committee approval was not needed for the study. The 
data we used for this study were publicly available aggregated 
data on annual numbers of THRs performed in Sweden and 
official statistics on population growth and projections of 
population growth. Registration of patients in the Swedish 
Hip Arthroplasty Register is regulated by Swedish law (SFS 
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Table 1 . Overview of the evolution of the Swedish population from 
1 970 to 201 0 and the official forecast for 2020 and 2030 



Life expectancy 



Year 


Population 


Foreign-born 


Men 


Women 


1970 


8,081,229 




72 


77 


1980 


8,317,937 




73 


79 


1990 


8,590,630 




75 


80 


2000 


8,882,792 


1,003,798 


77 


82 


2010 


9,415,570 


1,384,929 


80 


84 


2020 


10,200,459 


1,889,364 


81 


85 


2030 


10,660,644 


2,007,019 


83 


86 






Age distribution 






<40y 


40-59 y 


60-79 y 


a80y 


1970 


4,460,168 


2,028,518 


1 ,402,289 


190,167 


1980 


4,610,136 


1 ,866,990 


1 ,577,470 


263,341 


1990 


4,490,257 


2,150,431 


1 ,580,364 


369,578 


2000 


4,508,136 


2,401,292 


1 ,520,802 


452,562 


2010 


4,613,254 


2,450,925 


1 ,854,487 


496,904 


2020 


4,966,091 


2,570,152 


2,117,308 


546,908 


2030 


5,034,633 


2,595,266 


2,218,356 


812,089 



2008:355). Patients are informed about the registration by the 
hospital, no informed consent is needed, and patients may 
request withdrawal at any time. 

Data 

We obtained data for the period 1968-2012. We collected the 
annual number of THRs in Sweden from the Swedish Hip 
Arthroplasty Register (www.shpr.se). Annual numbers of 
THRs from 1968 to 1978 were based on data gathered from 
hospital records at the start of the registry. From 1979 to 1991, 
aggregated data on THRs were reported to the registry at the 
hospital level. Since 1992, all the relevant hospitals in Sweden 
have been reporting patient-specific data on all THRs to the 
registry with 98-99% completeness (Garellick et al. 2012). 
Publicly available data concerning the evolution of the Swed- 
ish population were downloaded from the Statistics Sweden 
website (www.scb.se). We downloaded historical data from 
1968 to 2012 and projections from 2013 to 2110— the total 
population, and the number of residents aged 40 years or 
older. From 1970 to 2010, the Swedish population increased 
by almost 1.5 million, there was a 7-year gain in life expec- 
tancy, and an increased proportion of elderly people (Table 
1). Forecasts by Statistics Sweden (from 2012) suggested that 
these trends would persist and that the proportion of elderly 
citizens would increase further in the coming decades. 

Statistics 

In order to adjust for differences in age distribution over the 
years, the annual incidence of THR was calculated for 10 5 
Swedish residents aged 40 years or older. This age group 
accounts for 99% of all THRs in Sweden. The estimated inci- 
dence served as outcome for the regression modeling while 



calendar year served as input. The regression models were 
built in order to forecast the incidence of THR operations per 
10 5 Swedish residents aged 40 years or older in the decades 
after 2012 and to estimate the maximum incidence per 10 5 
Swedish residents aged 40 years, if the model allowed. 

Two types of regression analysis were considered. First, we 
used Poisson regression analysis (Kurtz et al. 2007a, Patterson 
et al. 2009, Bini et al. 201 1). Poisson regression estimates the 
expected number of THRs per year and assumes a continuous 
growth. Thus, at least theoretically, the incidence can reach 
10 5 of 10 5 persons and if the results are used for projections, it 
could indicate unreasonably high numbers. 

Secondly, we adapted a regression framework that assumes 
the existence of an upper threshold, i.e. an asymptote that 
depicts the forecasted maximum incidence. The asymptote 
was estimated empirically from 3 competing models: asymp- 
totic, logistic, and Gompertz regression (Turner et al. 1969, 
Park and Lim 1985). Parameters of the Poisson regression 
were estimated via maximum likelihood. We used non-lin- 
ear least squares with the Levenberg-Marquardt algorithm 
to estimate the asymptote and parameters of the asymptotic, 
logistic, and Gompertz regressions (More 1978). 95% con- 
fidence intervals (CIs) were calculated based on the profile 
likelihood. 95% prediction intervals (Pis) were calculated 
based on error propagation using approximation by first- 
order Taylor expansion. 

Competing models were compared with the Akaike informa- 
tion criterion (AIC) and Akaike weights (wAIC) (Bozdogan 
1987, Wagenmakers and Farrell 2004). AIC penalizes model 
complexity and indicates which model best fits the data. As it 
lacks a direct and easy-to-interpret scale, AIC values are often 
transformed to wAIC. A model's wAIC is interpreted as the 
posterior probability of being the best model given the data 
and the set of tested competing models. Akaike weights were 
also used to calculate a weighted average of the estimated 
asymptote, thus obtaining estimates that incorporate model 
uncertainty (Lukacs et al. 2010, Symonds and Moussalli 
2011). Influential observations, whose removal would cause 
large change in the estimated asymptote, were detected with a 
jack-knife resampling. An observation was judged as influen- 
tial if it induced a change in the estimated asymptote by twice 
its standard error divided by the square root of the sample size. 

Statistical analyses were implemented in R 3.0.1 (Elzhov et 
al. 2013, R Core Team 2013). 



Results 

Between 1968 and 2012, 387,674 THRs were performed in 
Sweden according to the Swedish Hip Arthroplasty Register. 
Incidence in Swedish citizens aged 40 years and older and the 
number of THRs showed an increasing trend (Figure 1). While 
in 1968, only 179 THR operations were performed (with an 
incidence of 5 per 10 5 Swedish citizens aged 40 years and 
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Figure 1 . Annual number of total hip replacements registered in 
Sweden, 1968-2012. 



older), by 1970 this had increased to 906 (with an incidence of 
25). From 1970 to 1980, the increase was more than 5-fold — 
to 5,162 per year (incidence 139). Between 1980 and 1990, 
the annual number of operations almost doubled to 9,254 
(incidence 225). Between 1990 and 2000, the increase was 
lower and 1 1,329 operations were recorded in 2000 (incidence 
259). In 2010, 15,945 THRs were performed, giving an inci- 



dence of 332 per 10 5 Swedish citizens aged 40 years and older. 

Poisson regression offered a poor fit to the data, with a 
scaled Pearson chi-square statistic of 14.5 (1 represents a per- 
fect fit) and significant difference between the expected and 
observed data points (p<0.0001). Poisson regression forecast 
that the incidence of THRs will increase exponentially over 
the coming years, with a predicted incidence of 784 total hip 
replacements per 10 5 Swedish residents in 2030 and 1,133 in 
2040 (Figure 2A). 

The asymptote estimates from the competing asymptotic 
models varied considerably. The asymptotic regression was 
the most liberal (396, CF 256^460), the logistic regression 
was the most conservative (307, CF 289-332), and the Gom- 
pertz regression was in-between (324, CF 304-353). Model 
selection revealed the superiority of the asymptotic regression, 
with a posterior probability of being the best model — given 
the data and competing models — of 0.999. 

Asymptotic regression suggested a maximum incidence of 
THR of 396 per 10 5 Swedish residents aged 40 years and older 
(Figure 2B). This upper limit is forecast to be reached around 
the year 2107. Increasing incidence and population together 
with the increasing proportion of Swedish citizens aged 40 
years and older will lead to an increase in THR performed in 
the coming 2 decades (Table 2). 

Model diagnostics identified 4 influential points that would 
affect the estimated asymptote. The years 1987, 1991, 2009, 
and 2010 showed a sudden increase in THRs compared to pre- 
vious years. Removal of 1987 led to an estimated asymptote 
of 406; while removal of 1991 led to an estimated asymptote 
of 404. Removal of 2009 and 2010 both led to an estimated 



Incidence of Hip Replacement 




1970 1990 2010 2030 1970 1990 2010 2030 



A Year B Year 

Figure 2. The recorded and projected incidence of total hip arthroplasty per 10 5 Swedish residents aged 40 years or 
older, 2013-2040. A. Projections based on Poisson regression. B. Projections based on asymptotic modeling, with the 
gray horizontal line representing the estimated asymptote and the associated 95% confidence intervals of the best model; 
asymptotic regression. 
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Table 2. Prognosis for the evolution of the Swedish population between 2013 and 2030: the expected 
number and proportion of Swedish residents aged 40 years and older together with the predicted 
incidence per 10 5 Swedish residents aged 40 years and older and the expected number of total hip 
replacements (THRs) in Sweden as predicted by asymptotic regression 





Total 


Pnni ilatinn 

r uuuiqliui i 


Prnnnrtinn 

1 1 UUUI UKJ\ 1 






pYnpptpH nn 

1 A|JCTjLCU 1 \\J. 


Year 


noni Motion 


> 40 y 


> 40 y 


InriHpnrp 


95% PI 3 


of THR<; 

ul 1 1 1 1 lo 


2013 


Q rcjo 709 


4,947,192 


0.513 


324 


288-357 


16 021 


2014 


9 737 738 




0.513 


324 




16 318 


2015 


9 821 281 


5 042 1 1 8 


0.513 


329 


293-363 


16 595 


2016 


9,905,548 


5,082,444 


0.513 


334 


295-366 


16^854 


2017 


9,986,306 


5,120,677 


0.513 


334 


297-368 


17,104 


2018 


10,063,638 


5,156,449 


0.512 


336 


299-371 


17,343 


2019 


10,135,790 


5,194,795 


0.513 


339 


301-373 


17,588 


2020 


10,200,459 


5,234,368 


0.513 


341 


303-376 


17,834 


2021 


10,259,221 


5,269,706 


0.514 


343 


305-378 


18,063 


2022 


10,314,592 


5,304,478 


0.514 


345 


306-380 


18,288 


2023 


10,368,078 


5,337,914 


0.515 


347 


308-383 


18,505 


2024 


10,418,813 


5,372,586 


0.516 


349 


309-385 


18,725 


2025 


10,466,388 


5,410,267 


0.517 


350 


311-387 


18,952 


2026 


10,511,030 


5,448,911 


0.518 


352 


312-389 


19,180 


2027 


10,552,673 


5,487,684 


0.520 


354 


313-390 


19,407 


2028 


10,591,303 


5,531,168 


0.522 


355 


315-393 


19,648 


2029 


10,627,078 


5,575,927 


0.525 


357 


316-395 


19,892 


2030 


10,660,344 


5,625,711 


0.528 


358 


317-396 


20,152 



prediction interval. 



asymptote of 388. If all 4 outliers were removed at the same 
time, the estimated asymptote was 398. This could lead to a 
biased estimation of the asymptote; however, the values with- 
out the outliers were close to the asymptote value of the full 
model and well within the confidence interval. In addition, 
correction with jack-knife resampling reduced its value only 
by 1 patient, to 395. As there was unequivocal support for 
one model by AIC, model averaging changed the estimated 
asymptote only at the decimal level. 



Discussion 

Based on historical data, we attempted to forecast the future 
need for total hip replacement in Sweden. We found that the 
incidence of THR in Swedish residents aged 40 years and 
older will continue to increase until a projected upper inci- 
dence of around 400 total hip replacements per 10 5 Swed- 
ish residents aged 40 years and older will be reached at the 
beginning of the next century. This, in combination with the 
increase in and ageing of the Swedish population, will lead to 
a higher demand for this intervention and an increased burden 
on healthcare resources and hospital budgets. In the year 2020, 
the estimated incidence of THR will be 341 and in 2030 it 
will be 358. Using official forecasted population growth data, 
about 18,000 operations can be expected to be performed in 
2020 and 20,000 can be expected to be performed in 2030. 

Our findings are in line with similar studies showing a con- 
tinuous increase in joint replacements in past decades, and 
forecasting a continued increase in demand in the future. Kurtz 



et al (2005) observed a 46% increase in THRs in the USA 
between 1990 and 2002. This increase is projected to be 174% 
in 2030, imposing a substantial economic burden on health- 
care systems (Kurtz et al. 2007a, b). The historical increase 
in THR is not necessarily the result of a higher incidence of 
osteoarthritis, but merely an increase in operation rate (Skytta 
et al. 2011). Successful and widespread surgical treatment of 
osteoarthritis has a history of less than half a century, and non- 
surgical treatment is presumably still the most common choice 
globally. The epidemiology of conditions of the hip that may 
be treated with hip replacement varies in different populations 
(Makela et al. 2010). Hereditary and lifestyle factors as well 
as social, cultural, and economic factors will dictate how the 
demand for joint replacement surgery will change over time 
in different countries and populations. Accurate projections of 
future demand must consider that not every member of society 
needs or is eligible for THR. Modeling approaches of con- 
temporary projection studies have not considered this upper 
limit. Poisson regression, the workhorse of projection stud- 
ies, allows the projections to increase to infinity, at least theo- 
retically. This could lead to a serious overestimation of future 
incidence of total joint replacement. 

The methodological approach we adopted in the present 
work assumes that the incidence has an upper limit. Bini et 
al. (2011) reported a slowing demand for total joint replace- 
ment, with growth rates decreasing from 18% in 2002 to 3% 
in 2009. We observed a similar trend in the Swedish popu- 
lation. While the growth rate of THR was rather steep from 
the 1960s to the 1990s, after 2000 the growth slowed down. 
Barring the unexpected, such as medical breakthroughs, we 
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expect that the incidence of THR will continue to increase, 
even though the growth rate will be slowing continuously. In 
Sweden, primary osteoarthritis stands for around four-fifths of 
the diagnoses leading to THR. Thus, in the coming decades 
the incidence of THR will probably be determined to a large 
degree by the incidence of osteoarthritis. 

Historical data can help to forecast the trends in growth rates 
and to assess an upper level for the incidence of total joint 
replacement. The regression framework we considered allows 
estimation of an asymptote, an upper limit, and the associated 
confidence intervals. Estimation of the asymptote not only has 
a clinical/economic value but it is also important in obtain- 
ing a correct estimate of the growth rate. The drawback of 
the present framework is the increased computational com- 
plexity. Sample size requirements are likely to be higher than 
for the Poisson regression, and adjustment for other covari- 
ates is not straightforward. Calendar year used as predictor 
serves as a proxy for socio-economic characteristics of society 
in any particular year but cannot account for patient-related 
characteristics. Changes in the distribution of ages cannot be 
accounted for in these models, and this is a serious limitation. 
Nevertheless, we believe that the asymptotic framework that 
we adopted offers a biologically acceptable model and that 
its applicability in forecasting joint replacement is more plau- 
sible than Poisson regression. 

Naturally, every extrapolation from historical data — even to 
the immediate and near future — bears uncertainties (Davies et 
al. 2013). Independently of the statistical framework applied, 
the projections are based on trends inferred from historical 
data. These trends might change, and such trend changes are 
already palpable. Today, both younger, healthier patients and 
older, more frail patients undergo hip replacement, patients 
who would probably not have qualified for such an operation 
2-3 decades ago. An additional source of error is the change 
that the Swedish population goes through. At the end of 2012, 
Sweden had 9.6 million inhabitants, 14% of which were for- 
eign-born immigrants. This proportion was lower in previous 
decades (e.g. 9% in 1990), and it is likely to increase in the 
coming ones (the official forecast for 2020 is 22%, and for 
2030 it is 23%). Once operated, immigrants benefit from the 
operation as much as Swedish nationals (Krupic et al 2013). 
However, the extent to which immigrants take advantage of 
THR operations is not known. If the rate of THR varies among 
ethnic groups with the increasing heterogenization of the 
Swedish population, the validity of forecasting will decrease. 
Changes in lifestyle and longer life expectancy may increase 
the incidence of osteoarthritis, thus increasing demand for 
total replacement. Economic down-turns might constrain hos- 
pital budgets and reduce operation rates. Moreover, there may 
be new techniques or new medications to treat osteoarthritis, 
making total replacement obsolete. 

In conclusion, we expect that the incidence of total hip 
replacement in Swedish residents aged 40 years or older will 
slowly increase in the coming decades. The major drive behind 



the increased total number of operations will probably be the 
population increase and the growing proportion of residents 
aged 40 or more. 



SN participated in the design of the study, prepared databases, carried out 
the statistical analyses, and was responsible for drafting the manuscript. MG 
revised the manuscript. CR participated in the design of the study and revised 
the manuscript. OR supervised the study, participated in its design, and par- 
ticipated in drafting the manuscript. All the authors read and approved the 
final manuscript. 



The Swedish Hip Arthroplasty Register is funded by the Swedish Associa- 
tion of Local Authorities and Regions and by Region Vastra Gbtaland. These 
sources of funding had no role in study design, data collection and analysis, 
decision to publish, or preparation of the manuscript. No commercial funding 
was received for this study. 



Berger A, Bozic K, Stacey B, Edelsberg J, Sadosky A, Oster G. Patterns of 
pharmacotherapy and health care utilization and costs prior to total hip 
or total knee replacement in patients with osteoarthritis. Arthritis Rheum 
2011; 63 (8): 2268-75. 

Bini S A, Sidney S, Sorel M. Slowing demand for total joint arthroplasty in a 
population of 3.2 million. J Arthroplasty 201 1 ; 26 (6): 124-8. 

Bozdogan H. Model selection and Akaike Information Criterion (AIC)-The 
general-theory and its analytical extensions. Psychometrika 1987 52(3): 
345-70. 

Davies C, Briggs A, Lorgelly P, Garellick G, Malchau H. The "hazards" of 
extrapolating survival curves. Med Decis Making 2013; 33 (3): 369-80. 

Elzhov T V, Mullen K M, Spiess A N, Bolker B. minpack.lm: R interface 
to the Levenberg-Marquardt nonlinear least-squares algorithm found in 
MINPACK, plus support for bounds. R package version 1 1-7. 2013; http:// 
CRAN.R-project.org/package=minpack.lm. 

Garellick G, Karrholm J, Rogmark C, Rolfson O, Herberts P. Swedish hip 
arthroplasty register. Annual report 201 1. 2012; ISBN 978-91-980507-1-4. 

Kotlarz H, Gunnarsson C L, Fang H, Rizzo J A. Insurer and out-of-pocket 
costs of osteoarthritis in the US evidence from national survey data. Arthri- 
tis Rheum 2009; 60 (12): 3546-53. 

Krupic F, Eisler T, Eliasson T, Garellick G, Gordon M, Karrholm J. No influ- 
ence of immigrant background on the outcome of total hip arthroplasty. 
Acta Orthop 2013, 84 (1): 18-24. 

Kurtz S, Mowat F, Ong K, Chan N, Lau E, Halpern M. Prevalence of primary 
and revision total hip and knee arthroplasty in the United States from 1990 
through 2002. J Bone Joint Surg (Am) 2005; 87 (7): 1487-97. 

Kurtz S, Ong K, Lau E, Mowat F, Halpern M. Projections of primary and 
revision hip and knee arthroplasty in the United States from 2005 to 2030. 
J Bone Joint Surg (Am) 2007a; 89: 780-5. 

Kurtz S M, Ong K L, Schmier J, Mowat F, Saleh K, Dybvik E, et al. Future 
clinical and economic impact of revision total hip and knee arthroplasty. J 
Bone Joint Surg (Am) 2007b; 89: 144-51. 

Lukacs P M, Burnham K R Anderson D R. Model selection bias and Freed- 
man's paradox. Ann Inst Stat Math. 2010; 62 (1): 1 17-25. 

Makela K T, Peltola M, Hakkinen U, Remes V. Geographical variation in 
incidence of primary total hip arthroplasty: a population-based analysis of 
34,642 replacements. Arch Orthop Trauma Surg 2010; 130 (5): 633-9. 

More J. The Levenberg-Marquardt algorithm: Implementation and theory. In: 
Numerical Analysis. Lecture Notes in Mathematics (ed. G A Watson). 630: 
Springer Berlin Heidelberg; 1978: 105-16. 



Acta Orthopaedica 2014; 85 (3): 238-243 



243 



Park E W, Lim S M. Empirical estimation of the asymptotes of disease prog- 
ress curves and the use of the Richards generalized rate parameters for 
describing disease progress. Phytopathology 1985; 75 (1): 786-91. 

Patterson C C, Dahlquist G G, Gyiiriis E, Green A, Soltesz G. Incidence 
trends for childhood type 1 diabetes in Europe during 1989-2003 and pre- 
dicted new cases 2005-20: a multicentre prospective registration study. 
The Lancet 2009; 373 (9680): 2027-33. 

R Core Team. R: A Language and environment for statistical computing. R 
Foundation for Statistical Computing 2013; Vienna, Austria(http://www.R- 
project.org/). 

Skytta E T, Jarkko L, Antti E, Huhtala H, Ville R. Increasing incidence of hip 
arthroplasty for primary osteoarthritis in 30-to 59-year-old patients. Acta 
Orthop2011; 82(1): 1-5. 



Symonds M R E, Moussalli A. A brief guide to model selection, multimodel 
inference and model averaging in behavioural ecology using Akaike's 
information criterion. Behav Ecol Sociobiol 2011; 65 (1): 13-21 

Turner M E, Monroe R J, Lucas H L, Jr. Generalized asymptotic regression 
and non-linear path analysis. Biometrics 1961; 17 (1): 120-43. 

Turner M E, Jr., Blumenstein B A, Sebaugh J L. 265 Note: A generalization of 
the logistic law of growth. Biometrics 1969; 25 (3): 577-80. 

Wagenmakers E J, Farrell S. AIC model selection using Akaike weights. Psy- 
chon Bull Rev 2004; 11 (1): 192-6. 

Woolf A D, Pfleger B. Burden of major musculoskeletal conditions. Bull 
World Health Organ 2003; 81 (9): 646-56. 



